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Fermions play an essential role in many areas of quantum physics and it is desirable to understand 
the nature of entanglement within systems that consists of fermions. Whereas the issue of sepa- 
rability for bipartite fermions has extensively been studied in the present literature, this paper is 
concerned with maximally entangled fermions. A complete characterization of maximally entangled 
quasifree (Gaussian) fermionic states is given in terms of the covariance matrix. This result can be 
seen as a step towards distillation protocols for maximally entangled fermions. 
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I. INTRODUCTION 

Since fermions play an essential role in almost all ar- 
eas of quantum physics, it is natural to study the nature 
of entanglement within fermionic systems. The notion 
of entanglement relies on the bipartite or even multipar- 
tite structure of the underlying quantum system. We 
focus here on the bipartite case by considering two par- 
ties, called Alice and Bob. To determine a bipartite sys- 
tem, one has to specify the observables that are accessible 
within Alice's and Bob's laboratories. The collection of 
all possible measurements that can be made is described 
by the global observable algebra. Following the reason- 
ing of Refii, a bipartite system can mathematically be 
described by a pair of local observable subalgebras sitting 
inside the global observable algebra. The two local alge- 
bras represent the observations that can be made within 
Alice's and Bob's laboratories, respectively. The essen- 
tial requirement is that each measurement performed by 
Bob cannot disturb the measurements made within Al- 
ice's laboratory. Therefore, each of Alice's operators has 
to commute with all of Bob's operators. Once a bipartite 
system is fixed by choosing the local algebras, the no- 
tion of entanglement can be introduced relatively to this 
choice: product states are characterized by the property 
that the expectation values of a product of one of Alice's 
operators and one of Bob's operators factorize. Taking 
the convex hull of all product states yields the set of sep- 
arable states and the issue of entanglement is to detect 
those states which are not separable. We refer the reader 
to Refsj^i^, where the characterization of separable states 
on fermion systems is discussed. 

In comparison to the considerations in Refi^, we fo- 
cus here on the "converse" of separability, namely the 
characterization of maximally entangled fermions. When 
dealing with fermions, the main issue is that, due to 
the canonical aniticommutation relations (CAR) of Fermi 
fields, generic operators of Alice and Bob do not commute 
among each other. As we shall see, the standard bipar- 



tite structure discussed above can be restored either by 
restricting the analysis to even products of Fermi fields, 
as mainly done in this paper, or by introducing a twisted 
tensor product structure which automatically accounts 
for the alternating signs related to the commutation of 
Fermi fields belonging to different subsystems. 

The restriction to even products of Fermi fields corre- 
sponds to the requirement that physical observables pre- 
serve the parity of the particle number. Since correlation 
experiments for the detection of entanglement are built 
from observables that belong to the local subsystems of 
Alice and Bob, it is natural to restrict the analysis to 
operators which preserve not only the global parity, but 
also the local parity. In mathematical terms, we will 
then exploit the restriction to the even-even part of the 
global algebra. Consequently, we shall speak of accessible 
entanglement, relative to the restriction under consider- 
ation. Clearly, once the even-even part of the algebra is 
considered, operators belonging to different subsystems 
commute among each other, giving rise to a tensor prod- 
uct structure analog to the qubit case. 

As far as this restriction is concerned, our main result 
is the full characterization of quasifree (i.e., Gaussian) 
maximally entangled states. We prove that, in order to 
be maximally entangled, a fermionic quasifree state must 
have a covariance matrix of a well-defined form, which is 
unique up to local Bogolubov transformations. 

This analysis obviously requires a precise definition of 
fermionic maximally entangled states in the presence of 
the above restriction. To this purpose, we exploit the 
conservation of local parity to further decompose the op- 
erator algebra. Indeed, it turns out that the structure 
of Alice's and Bob's observable algebras corresponds to 
the following situation: Alice possesses a direct sum of 
matrix algebras of the form 0^ 'Qi^^'^x whereas Bob pos- 
sesses the algebra 0^ 1x®'^x, where Sl^, = = Md{C) 
are full matrix algebras. Here, the label x refers to the 
number parity of the two subspaces (four combinations 
are possible, i.e., H — h, H — , — +, , where + and — cor- 
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respond to even and odd. respectively), and, assuming 
that both Ahee and Bob own n fermionic modes, one has 
d ~ 2"^^. The direct sum structure of the observable al- 
gebras can be interpreted as a super- selection rule, where 
the label x plays the role of a conserved charge or quan- 
tum number. Therefore, we are faced here with the issue 
of entanglement in the presence of a super-selection rule. 

While correlation experiments for the detection of en- 
tanglement only involve observables belonging to the 
quasi-local observable algebra given by 0^ Sl^, ® 05^ , the 
fermionic state g under investigation can be prepared in 
a more general way and can be thought of as belonging 
to the global observable algebra. ^ The restriction of g to 
the even-even part of the algebra, i.e., to the quasi- local 
observable algebra, is given by a density operator g'^'^ 
of block-diagonal form, p++ ~ where is a 

normalized probability distribution and gx is a density 
matrix in 21^: <Ei ^x- Pure states are given by rank-one 
projections ax in ^x ® ^x and the situation is effectively 
the same as for a simple tensor product since there is 
only one non-zero term that contributes to the direct 
sum. A suitable entanglement measure is the entropy 
of entanglement E{ax) = S{tT:<s^{'^x)) where iT:<B^{'Jx) is 
the reduced density matrix on Alice's subblock St^ and 
S{p) = — tr[plog2(p)] is the von Neumann entropy. If we 
consider a density operator of the form g^'^ = ®xP^'^^ 
with rank-one projections Ux, then (due to the supcrsc- 
lection rule) this is the unique convex decomposition of 
g"*""*" into pure states. For this class of states the entangle- 
ment of formation Ep can easily be computed. It is just 
given by the mean value of the entropy of entanglement 
of each block, i.e., Epig'^^) = ^xP^^i'^^)- we can 
ask for those states which maximize the entanglement 
of formation. They are given by all density operators 
of the form g"*""*" = ^^Pxt^x where (Tx is a projection 
onto a maximally entangled vector. For all these states 
we find for the entanglement of formation the maximal 
value Ep{g'^~^) = log2(c?). This is precisely the notion of 
maximal entanglement that we apply here to the inves- 
tigation of entanglement in fermionic systems. 

The discussion of fermionic entanglement can also be 
applied to investigate entanglement in spin chain sys- 
tems with the help of the Jordan- Wigncr transformation. 
Then, preservation of global parity corresponds to preser- 
vation of total magnetization. Since magnetization does 
not give rise to superselection rules, a fermionic system 
formally arising from the Jordan- Wigner transformation 
of a spin chain is not necessarily subject to the physi- 
cal constraints discussed above. In particular, the two 
local observable algebras arc not restricted to the even 
part, given by even products of Fermi fields, and their 
commutativity is not guaranteed. In such a context, it 
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^ Differently from the quasi-local observable algebra, which pre- 
serves local parity, the global observable algebra only preserves 
global parity and therefore includes operators which cannot be 
expressed as products of local observables. 



FIG. 1: The figure illustrates a spin chain that consists of two 
blocks of spins. Alice controls the block {1, 2, 3, • • • ,7} (dark 
grey squares) and Bob possesses the block {8, 9, ■ ■ ■ , 14} (light 
grey squares). The black lines are connecting those positions 
which are prepared in a maximally entangled two qubit state. 



is useful to consider the twisted tensor product structure 
mentioned above, which takes into account the anticom- 
mutation relations of Fermi fields similarly to the Jordan- 
Wigner transformation. 

The outline of the paper is the following. We start by 
giving an elementary overview (Section [Tl| of the basic 
ideas, and discuss, as an instructive example, the sim- 
plest non-trivial case which corresponds to two fcrmion 
modes per part. We state our main result on the char- 
acterization of quasifrce maximally entangled fcrmions 
in terms of this example first. The structure of general 
bipartite fcrmion systems in finite (but arbitrary large) 
dimensions is presented in Section [1111 We take here ad- 
vantage of Araki's self dual formalisn:^^, which enables 
to write many formulas in a more handy form. The main 
result on the characterization of maximally entangled 
quasifree fermion states is derived in Section IIVI We 
give here an explicit formula for the covariance matrix 
for a maximally entangled quasifree fermion state. We 
also show that, up to local Bogolubov transformations, 
there is a unique maximally entangled quasifree fermion 
state. Bogolubov transformations are reversible physical 
operations that arc characterized by the property to map 
quasifree states into quasifree states (see Refsi^i^iii^ for 
general quasifree completely positive maps) . We also dis- 
cuss a one-parameter group of (global) Bogolubov trans- 
formations (which can be interpreted as the unitary time 
evolution due to a proper interaction) mapping a product 
state into a maximally entangled fermion state for a suit- 
able value of the associated parameter (time) . As a fur- 
ther example, we consider in Section [Via- spin chain sys- 
tem consisting of two equally sized blocks of spins. The 
system is prepared in a state where each site of one block 
is maximally entangled with a companion site of the other 
block, according to Fig. HI Applying the Jordan- Wigner 
transformation the corresponding state on the fermion al- 
gebra is quasifree and maximally entangled in the sense 
of our basic definition. Finally, in Sec. IVII we provide 
the description of the twisted tensor product structure 
which is necessary to extend the discussion beyond the 
even-even subalgebra. 
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II. OVERVIEW 

In order to have a closer look at the structure of Al- 
ice's and Bob's observable algebras, we consider first the 
most simple example which consisits of two fermionic 
modes per part. The creation and annihilation opera- 
tors ct , d (generically denoted by cf in the following) 
are given here in the "occupation number representa- 
tion". By fixing an ordering of the modes first, this rep- 
resentation is given on the 2^ = 16-dimcnsional Hilbcrt 
space Ti. that is spanned by basis vectors \q1q2q3q4:) ~ 
(4)«i(4)'?'(4)'^'(ci)'"|0000), labeled by binary strings 
of length four, i.e., gi = 0, 1. The 16-dimcnsional Hilbcrt 
space Ti. can be decomposed into a direct sum of four 
4-dimensional subspaces. Each vector E H can be rep- 
resented by a block decomposition 



(1) 



where the four block vectors "if^y G 



with 



x,y — ±, arc simultaneous eigenvectors of the parity 
operators of both subsystems. For example, ^'++ G 
n++ = span{|0000), |1100), 10011), |1111)}. Alice's ob- 
servable algebra is generated by the operators cfc^" with 
I, j = 1, 2. If we apply these generators to the basis vec- 
tors, it is easy to see that Alice's observables respect the 
block decomposition Eq. ^ by acting trivially on Bob's 
subspace and preserving the parity of Alice's subspace. 
For example, cfcf^++ G H++ for ij = 1,2. 

Putting all these together, a general operator of Alice's 
local observable algebra can be represented by the block 
matrix 



A = 



A_ (g) 1. 



(2) 



with 2x2 matrices A± G Af2(C). By an analogous rea- 
soning, a general operator that belongs to Bob's local 
observable algebra has the form 



B = 



\ 



1+ ® B. 



1- ® B4 



\ 



1_ «) B_ / 



(3) 



with 2x2 matrices B± G Af2(C). This representation of 
the local observable algebras shows directly the fact that 
Alice's operators commute with Bob's. 

We suggest here a notion for maximally entangled bi- 
partite fermion states which is related to the maximiza- 
tion of the entanglement of formation. In view of our 
example above, we declare a state 



Q = 



(4) 



where Exy, x,y = ±, are projections onto maximally 
entangled vectors VL^y, to be maximally entangled. 

The example given above, is concerned with two 
fermion modes per site (two for Alice and two for Bob) 
and the total dimension of the Hilbcrt space is 2^. More 
generally, if each of Alice and Bob possesses n fermion 
modes, then the dimension is 2^". Therefore the dimen- 
sion of the Hilbcrt space of the system increases exponen- 
tially with the number of fermion modes. It is therefore 
highly desirable to find ways of treating at least some as- 
pects of large fermion systems without actually having to 
write out state vectors component by component. Here 
the class of quasifree fermion states is a suitable choice. 
They can be characterized by simple combinatorial data, 
which do not grow exponentially but at most quadratic 
with the number of modes. Quasifree states are deter- 
mined by the two-point correlation fmictions. i.e., the ex- 
pectation values (cfc^). For 2n fermion mode a quasifree 
state is determined by a 4n x An matrix which have to 
fulfill suitable constraints that are related to the positiv- 
ity of the density operator for the state. Although the 
mathematical structure of quasifree states can be treated 
easily, this class of states is sufficiently complex to tackle 
many non-trivial physical systems, in particular to find 
ground states of interacting spin chains with the help of 
the Jordan- Wigner transformation. 

To determine a quasifree state, we arrange the expec- 
tation values (cfcf), i,j = l,...2ri, into the 4n x 4n 
hermitian matrix 



S 



Skk Sub 
Sbk Sbb 



(5) 



where Sba = SIb- Any of the 2n x 2n submatrices 5'xy, 
with X, Y = A, B, has the structure 



Sxv — 



c-- c — h 
"-■XY '^XY 



and the n x n submatrices S'xy', with x,y = 



(6) 



are 



(cjcfe), {SIb )jfc = (c]c„+fe), and so 



given by (S'+ )jk 

on. Here, (cfcf) ; ..^ 
for the case n = 2, we fully list the off-diagonal block S^s,^ 
which plays the essential role in view of entanglement: 



ii^iQscfcf). For the sake of clarity. 



Sab — 



( (4C3> 
(4c3) 
(C1C3) 

V (C2C3) 



(4C4) 
(C1C4) 
(C2C4) 



(44) 
(44) 
(ci4) 

(C24> 



(44) \ 
(44) 

(ci4> 

(C24> / 



(7) 



The matrix S is called the covariance matrix of the 
quasifree state and we denote by gs the corresponding 
density matrix. The diagonal blocks of S are related to 
the Alice-Alice and Bob-Bob correlations. The ordering 
of the entries of S is chosen in agreement with Araki's 
sclfdual formalism, which is explained in the next section. 

The main result of this paper gives a complete char- 
acterization of maximally entangled quasifree states for 
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the case that Alice and Bob possess an arbitrary num- 
ber n of fermion modes by giving an exphcit form of the 
covariancc matrix. For two modes per site, the standard 
example for a maximally entangled fermion state is given 
by the covariance operator 



2 I -il4 



il4 

l4 



(8) 



where here and in the following the symbol denotes 
the n X n identity matrix. It is not difficult to check that 
P is a projection which implies that the quasifree state 
gp is pure on the full fermion algebra that is generated by 
cf ,cf ,cf ,cf . The first thing one observes here is that 
the diagonal blocks Paa = Pbb = 1/2 are given by one- 
half times the unit operator. Thus the state restricted 
to Alice's fermion operators has covariance 1/2. As we 
shall prove, the quasifree state with covariance 1/2 is the 
totally mixed state, i.e., the corresponding density matrix 
is ~ 1/4- Recall that Alice's fermion operators 

cf,cf generate the algebra Af4(C) of all complex 4x4 
matrices. The fact that the reduced density matrix is 
maximally mixed is the essential ingredient to prove that 
the reduced density matrix ofgp to the subalgebra that 
is generated by Alice's and Bob's observables has indeed 
the form of Eq. 



III. BIPARTITE FERMION SYSTEMS AND 
THE CAR ALGEBRA 

Let us now consider the most general case, where Alice 
controls a finite set A of n fermion modes, whereas Bob 
possesses another set B of m fermion modes. We denote 
by the full local algebra of Alice, which consists of all 
possible linear combinations of products of creation and 
annihilation operators cj, Ca belonging to her set of modes 
a G A. Alice's local observable algebra is instead given by 
the even part of 5^a, which is denoted by and is given 
by linear combinations of products of an even number of 
Fermi operators. The same notation is used for Bob, by 
an obvious replacement of the subscripts. The full global 
algebra generated by Alice's and Bob's fermion operators 
together is denoted by S^ab- The global observable algebra 
is the even part of 3^ab and is denoted by 3^^. Finally, the 
even-even part of the global algebra, corresponding to the 
tensor product ^'l ^ dt ^ is denoted by dts^ ■ Similarly 
to the even part, one can also define the odd part of an 
algebra, e.g., tJa" : i?b' : i?ab • Note however that, differently 
from the even part, the odd part is only a linear subspace 
and not an algebra, as products of odd operators yield 
an even operator. 

The bipartite structure is given by the natural separa- 
tion of fermionic modes relative to the (full) subalgebras 
5'a and JJb- Fermionic entanglement is based on the prop- 
erties of correlation functions (AB) = tr{gAB), i.e., the 
expectation values of operator products AB of a state g, 
where A belongs to Alice's fermion algebra 3^a and B is an 



operator in Bob's fermion algebra 5^b- In comparison to 
the standard entanglement theory, one is faced with the 
problem that the operators A and B do not necessarily 
commute. Namely, they anticommute whenever both A 
and B are odd operators, i.e., A e and B G 5^^. This 
can be seen as a non-local behaviour which makes the 
fermionic bipartite structure different from the standard 
tensor product. The appropriate description of this fea- 
ture is provided by the twisted tensor product discussed 
in Sec. IVII On the other hand, in the physical situa- 
tions where local parity is preserved, it is sufficient to 
consider the restriction to the even-even part of the al- 
gebra, where operators belonging to different subsystems 
commute and the standard tensor product can be used. 
This will be our assumption in Sees. IIVI and fVl In the 
rest of this section we will introduce a convenient nota- 
tion to deal with fermionic fields and fermionic Gaussian 
states. 



A. Araki's selfdual formalism 

In order to give a clear presentation, wc make use of 
Araki's selfdual formalism which enables to write formu- 
las in a more compact manner. The basic idea is to com- 
bine the creation and annihilation operators into one field 
operator by the linear combinations 



B{f) 



jGAB 



-4 



(9) 



The 2{n + m) complex coefficients /j are put into one 
vector / of a Hilbert space /Cab = C" © C" © © 
which we call the Hilbert space of the global one-particle 
system. Each vector / G i^AB can be decomposed into 
four blocks / = {f^, f^, f^, f^) where ft G C" are the 
one-particle components for Alice and /g G C™ are the 
one-particle components for Bob. Recall that Alice pos- 
sesses n = |A| modes and Bob has m = |B| modes under 
control. We introduce a complex conjugation F which al- 
lows to write the canonical anticommutation relation in a 
more elegant form. Making use of our block decomposi- 
tion, the complex conjugation is an antiunitary operator 
that can be represented by a 4 x 4 block matrix 



(10) 



where Jk (in our case we have k = n, m) is the com- 
plex conjugation on C*^ that replaces each component 
within a vector by its complex conjugate. By using the 
anticommutation relations for the standard creation and 
annihilation operators cj , Cj allows to write the canonical 
anticommutation relation in the following compact form: 





Jn 








Jn 




















Jrn 


V 





J m 






{BU).B{g)}^{Tf,g)t 



(11) 
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where ( , ) denotes the usual complex scalar products be- 
tween vectors. Moreover, the adjoint of B{f) can simply 
be computed according to 



B(/)t = B{Tf) 



(12) 



Indeed, the global fcrmion algebra 5^ab is completely de- 
termined by requiring that the field operator B{f) is 
complex linear in the one-particle vector / and by the 
relations Eqs. ([TT]) (fT2)) . 



B. Quasifree states 

Since S^ab is a finite dimensional full matrix algebra, 
the states on the bipartite fcrmion system are given by 
density operators g assigning to an operator X the ex- 
pectation value {X) = tr{gX). 

A quasifree state g, also called Gaussian fermion 
state, is completely determined by the two-point 
correlation functions, i.e., the expectation values of 
tr(gCjCi), tr(gcjcj), tr(£iCjcJ), tr(pCjCi) Araki's selfdual 
approach is very handy by putting all this expectation 
values into a single big covariance matrix. Each quasifree 
state g = gs '^s in one-to-one correspondence with a co- 
variance matrix which is a linear operator 5* on Ki^s that 
is between and 1 (0 < S* < 1) with respect to the oper- 
ator ordering and which fulfils the additional constraint 
S + rSr = 1. The expectation values of the state u; 
are related to its covariance matrix S by the following 
condition on the two-point correlation function: 



tT{gsBif)Big))^(Tf,Sg) . 



(13) 



The entries of the matrix S arc related to pair correlation 
functions as described in Sec. IIIII All higher correlation 
functions can be expressed in terms of sums of products 
of two-point functions according to Wick theorem, where 
only the expectation values of an even product of Fermi 
field operators are non vanishing (i.e., gs € Sts)- It is 
well known that a quasifree state is pure if and only if its 
covariance matrix S* is a projection, called basis projec- 
tion. One simple example for a basis projection can be 
given in terms of the 4x4 block matrix 



E := 



/ 1„ 

1™ 
V 



(14) 



according to the four-block-decomposition of /^ab- In- 
deed, a straightforward computation shows that the con- 
dition E + VET = 1 is fulfilled. The expectation val- 
ues for the standard creation and annihilation operators 
are ti{gEclca) = trlgEclcb) = 1 for a G A and b E B. 
The expectation values of all other possible combina- 
tions of Fermi operators are vanishing. In particular, 
tT{gEB{Qf,f)B{Qsh)) ~ which implies that gs is a 
product state, i.e., if we consider a product AB of an 



operator A from Alice's local fermion algebra and an op- 
erator B from Bob's local fermion algebra, then the ex- 
pectation value tr{gEAB) = tT{gEA)tT{gEB) factorizes. 



IV. MAXIMALLY ENTANGLED STATES ON 
BIPARTITE FERMION SYSTEMS 

In the following wc assume that Alice and Bob con- 
trol the same number of fermion modes n = m. Due 
to the bipartite nature of the fermion algebra 3^ab there 
is a natural notion of product states. A state g on 
3^AB is called a product state if the expectation value 
tr{gAB) = tr{gi^A)ti{gBB) factorizes for any pair of op- 
erators A, B belonging to Alice's and Bob's local fermion 
algebras, respectively. ^ The quasifree state given by 
Eq. (fTi|) is a particular example. As for standard entan- 
glement theory, a state on J^ab is entangled if it can not be 
written as a convex combination of product states^i^. But 
what about maximally entangled states? There is a nat- 
ural point of view for discussing entanglement in fermion 
systems. If we think of the observables that can be mea- 
sured within Alice's and Bob's local laboratories, we have 
to restrict to the even parts of the local subalgebras, i.e., 
and . In particular, each correlation experiment 
for detecting the presence of entanglement can only be 
built from observables that belong to the local even parts. 
As a consequence, the accessible entanglement of a bipar- 
tite fermion state g is the entanglement of the restriction 
to the even-even subalgebra S^^"*" which can be directly 
identified with the tensor product ^ ■ 

As already mentioned, the algebra ^Jab can be identified 
with a full matrix algebra of all linear operators on a 
Hilbcrt space with dimension 2^" . We choose now an 
identification that is most compatible with the bipartite 
nature, namely we identify JJab = ^{'H <8) H) with the 
matrix algebra of all linear operators on H (E)'H, where 
H F_(C") is the antisymmetric Fock space over C". 
We identify the set A of labels for Alice's fcrmion modes 
with A = {oi, 02, • • • , a„} and the set of labels for Bob's 
fermion modes with B = {bi, ■ ■ ■ ,bn}- Let Cj, Cj denote 
the standard creation and annihilation operators on Ti.. 
Then Alice's fermion operators act on (g) by Caj ~ 
Cj (g) 1 and Bob's fermion operators act on Ti (g) Ti by 
Cbj = 0®Cj, where 6 ^ (—1)-'^ is the parity of the particle 

number TV = X]j=i cjcj on Ti. The Fock space Ti = 
Ti+ © Ti_ decomposes into the even and odd particle 
number subspaces Ti± that correspond to the cigenspaces 
of 9 with eigenvalues ±1. The projections onto Ti± are 
1± = {t±9)/2 and clearly 1 = 1+©1_. Thus the tensor 
product 'H®'H can be expanded into a direct sum of four 
subspaces Ti (g) Ti = (Ti+ ® Ti+) ® (Ti+ ® Ti_) ® (Ti_ ® 
Ti-)-)0 (Ti_ (K)Ti_). The local parity number operators for 



^ Possible variations of this definition of product state arc dis- 
cussed in Ref™. 
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Alice and Bob are given by 8a = 



1 and Bb = 1 8) 6*. 

The even-even subalgebra ^^b' consits of all operators 

X commutes with 



X that preserve the local parity, i.e. 
Oa and Ob 

diagonal X = X++ ® X+_ ® X_+ This yields 

the desired decomposition 



Therefore the operator X e ^'a'b^ is block 



"AB 



^(Ha)®^(Hb);. 



(15) 



a,b=± 



where 7i+ and H- have the same dimension 2"~^. 

If the global fermion system is prepared in a state 
g, then the restriction g'^^ to the even-even part has 
a unique convex decomposition 



= ^ Pab CTab 

a,b=± 



(16) 



where aab is a state on ^{Ha) <X) ^{Hb)- We suggest 
here to characterize maximally entangled fermion states 
by the following definition: 

Definiton 1 A bipartite fermion state g on S^ab is called 
maximally entangled if each state aab within the decom- 
position Eq. p6p is maximally entangled in the usual 
sense. 

This point of view is related to searching for those 
states for which the restriction g"*""*" to the even-even part 
maximizes the entanglement of formation Ep. Recall 
that it can be evaluated according to 



Epig++) 



inf 



(17) 



where the infimum is taken over all possible decomposi- 



tions g 



++ 



'^Ix'^x into pure states ax and E denotes 



the entropy of entanglement. If the state is maximally 
entangled according to Dcf. [U the states aab inside de- 
composition pop , being maximally entangled in the usual 
sense, must be pure. Then, the direct sum structure of 
S^AB^ implies that the decomposition ITBl) is the only de- 
composition into pure states and hence 



E, 



<e^^) = E{aab) = log2 2"-i =n-l. (18) 



a,b=± 



The density matrix Uab is the projection onto a maxi- 
mally entangled vector Vlab in Jia ® Ti-b- This naturally 
induces a vector f2 = (9a,b=±\/Pab^ab in the full Hilbert 
space Ti®T-l which corresponds to a pure state g = 
on the global fermion algebra ^{H^Ti.) = 5'ab- Restrict- 
ing the analysis to the accessible entanglement mentioned 
above, yields a definition that only constrains the form of 
. The global fermion state g is then not uniquely de- 
termined. It is however natural to require it to be pure, 
as for usual maximally entangled states. Therefore, we 
will assume each maximally entangled bipartite fermion 
state to be pure on the global fermion algebra. 



A. The structure of maximally entangled quasifree 
states 

Concerning quasifree states, we are faced with the fol- 
lowing problem: Is it possible to characterize explicitly 
the set of those covariance matrices which correspond to 
a maximally entangled fermion state? 

Theorem 2 A quasifree pure state gp on the global 
fermion algebra ^ab with covariance matrix P is maxi- 
mally entangled if and only if there exists a unitary oper- 
ator UnB that maps Bob's one-particle space onto Al- 
ice 's one-particle space and which anticommutes with 
the complex conjugation FaJ/ab = ^Ui.'^T-g, such that the 
covariance operator is given by the 2x2 block matrix 



P = 



1 



U C/ab 



(19) 



Proof: Suppose a maximally entangled quasifree fermion 
state gp with covariance matrix P is given. Since gp is 
pure, the covariance P is a basis projection and the cor- 
responding density operator pp is the projection onto a 
vector Hp ~ j^^j. y/pab^ab with maximally entangled 
vectors fiab- As a consequence, the restriction to Alice's 
even subalgebra is a trace, i.e., for each Ai,A2 € 
we have tr(gp^i^2) = tr(ppA2Ai). Moreover, the re- 
striction of the quasifree state to Alice local fermion al- 
gebra 5^A is again a quasifree state with covariance op- 
erator QaPQh, where Qa is the projection onto Alice's 
one-particle space. Recall that each quasi free state is 
determined by the two-point correlation function. Ac- 
cording to Theorem [HI which we prove in the appendix, 
the covariance operator for a quasifree state which is 
a trace on the even part must be 1/2. Thus we con- 
clude that QaPQa = Qa/2. The same also holds for 
the restriction to Bob's even subalgebra and we also find 
QbPQb = Qb/2. With respect to the Ahcc-Bob split 
Kab = Ka(BKb into the local one-particle spaces, the co- 
variance operator can be written as a 2 x 2 block matrix 



P 



Ia Uab 
Uba 1b 



(20) 



Here Uab is an operator from Bob's one-particle space Kb 
to Alice's one-particle space Ka and J7ba maps conversely 
A'b to A'a. Since P is sclfadjoint, the off-diagonal blocks 
must fulfill L^AB — Uba- Moreover, P is idcmpotent, i.e., 
P^ ^ P which implies 



P' = 



1 / 1a 



4 

1 ( 1a 

2 I Uba 



f UabUba 
2C/ba 

Uab 
1b 



2Uab 

+ UbaUab 



(21) 



Comparing the diagonal blocks yields the two identities 
UabUba = UabUIb = 1a and UbaUab = uIqUab = 1b which 
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shows that ^/ab is a unitary operator. Finally we have to 
take care of the constraint P + TPT = 1 for P being a 
valid covariance matrix. With respect to the Alice-Bob 
split of the one-particle space, the complex conjugation 
r is block diagonal 



r = 



Fa 
Tb 



(22) 



Recall that by comparing this 2x2 block representation 
with the 4x4 block representation Eq. pO|) the complex 
conjugation on Alice's (Bob's) one-particle space is given 
by the 2 X 2 block matrix 



Jn 

Jr.. 



(23) 



where J„ is the natural complex conjugation on C". Ex- 
pressed in terms of the Alice-Bob split, the condition 
P + rPr = 1 becomes 



rpv = -( ^^U,^T^ 
2 ^ TbUIt, 1b 



-ul 



(24) 



= 1 - F 



By comparing the off-diagonal blocks, we obtain the de- 
sired relation Uic^T^ = — FaJ/ab- 

It remains to prove the converse. Suppose a unitary 
t^AB is given such that C/abFa = — EbC/ab holds. Then 
Eq. defines a basis projection P that corresponds to 
a pure state Qp = |rip)(ilp| with ilp ^ Ti ®Ti. With 
respect to the direct sum decomposition into even and 
odd particle number subspaces 7i± the vector VL p can be 
written as flp = ® j,f, ^/Pab^ab where flab are normalized 
vectors in 'Ha(S)'Hb- Since a quasifree state vanishes on the 
odd part of the global fermion algebra, the vector flp has 

to be an eigenvector of 9. Thus either = p ^ =0 

or = p = 0. Furthermore, the restriction of gp 

to Alice's local fermion algebra 5^a is a trace since, by 
construction, QuPQn = Qk/2. Since = [^{H+) ® 
a] e [^(H+) (E) 1], the restriction of £)p to £§{H±) (g) 1 is 

a trace. If p^ = p y = 0, then the restriction of of gp 

to 

is supported on the blocks ^(7i+) (g) 1+ and 3§{Ti-) ® 
]1_. On the other hand, if p++ = p-- = 0, then the 
restriction of of gp is supported on the blocks ^{H+) (8) 
I- and S§{Ti-) ® 1+. Therefore, the restriction of the 
pure state Uab = \^ab){^ab\ is a trace on £§{Tla) <8) It 
which implies that flab is a maximally entangled vector 
in 'Ha®'Hb and is a maximally entangled fermion state 
in the sense of Definition [TJ ■ 

In summary, the theorem tells us that all maximally 
entangled fermion states can be characterized by a uni- 
tary operator J7ab from Bob's one-particle space onto Al- 
ice's onc-particle space that anticommute with the com- 
plex conjugation Fa^/ab = — C/abFb- Since, according to 



our assumptions, Alice's and Bob's one-particle spaces 
have the same dimension 2n they can be canonically iden- 
tified /\a = K-& = C^". The simplest choice for the uni- 
tary J7ab is just given by J7ab = i^2m where l2n is the 
unit operator on C^". The global one particle space is 
given by C^" ® C^" and the basis projection simply looks 
like 



P.i 



l2n il2n 
— il2n l2n 



(26) 



We call the corresponding quasifree fermion state the 
standard maximally entangled fermion state. 



B. Local Bogolubov transformations 

As we will see soon, each quasifree maximally entan- 
gled fermion state can be transformed into the standard 
maximally entangled fermion state by local Bogolubov 
transformations. A Bogolubov transformation is given 
by a unitary operator V on Xab that commutes with the 
complex conjugation = V^F. Indeed, if we trans- 
form the Fermi field operators by B{f) i-^ B{Vf) the 
anticommutation relations are preserved. Namely, we 
have {B{Vf),B{Vti)} = {TVf,Vh)t = {Vrf,Vh)t = 
{Tf,h)l. This implies that B{f) B{Vf) extends to 
a so called Bogolubov automorphism f3v of the fermion 
algebra S^ab- Since we are concerned with finite di- 
mensions, there exists a unitary operator Uy such that 
/3v{B{f)) = B{Vf) = UlB{f)Uv. Bogolubov auto- 
morphisms map quasifree states onto quasifree state. In- 
deed, a quasifree state gs with covariance S is mapped 
to the quasifree state UygsUy = gvsv^ with covariance 

vsv^.^. 

A local Bogolubov transformation is block-diagonal 
with respect to the Alice-Bob split of the one-particle 
space, i.e., there are unitary operators 14 on K/^ and Vb 
on such that 



V 



V^A 
Vn 



(27) 



where and Vb commute with the local complex conju- 
gations Fa and Fb, i.e., they are Bogolubov transforma- 
tions on the local one-particle spaces K/^ and K^, respec- 
tively. 

By construction, a local Bogolubov transformation 
preserves Alice's and Bob's even subalgebras. Namely, 
the block-diagonal structure of V implies that an opera- 
tor B{fii) of Alice's local fermion algebra is mapped to 
B{V fi^) = B{VkfA) which is also contained in Alice's lo- 
cal fermion algebra. The corresponding statement also 
holds for Bob. Thus we can restrict the Bogolubov auto- 
morphism f3v to an automorphism f3y~^ to the even-even 
subalgebra 5^^"*" which, as mentioned above, can be di- 
rectly identified with the tensor product (E) dt of ^i" 
ice's even part and Bob's even part. Then, the automor- 
phism f3y~^ splits into a tensor product (3y^ — (iy^ ® Py^, 
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where jSy^ and (3y-^ are the restrictions of the Bogolubov 
automorphisms and Pvb to Ahce's and Bob's even 
subalgebra. Hence a local Bogolubov transformation cor- 
responds to a local operation in the usual sense. 

If we apply a local Bogolubov transformation V to the 
standard maximally entangled state whose basis projec- 
tion is given by the Pst of Eq. ((26|) . we obtain a new 
maximally entangled state with basis projection VPstV^ 
which has the 2x2 block decomposition 



t 



l2n 



l2n 



(28) 



Let J/ab be a unitary that anticommutes with the com- 
plex conjugation and therefore determines a maximally 
entangled fermion state. Then — IL^ab commutes with the 
complex conjugation and 



W 



-i(7AB 

l2n 



(29) 



is a local Bogolubov transformation that acts on Bob's 
system alone. If we apply it to the standard maximally 
entangled state, then we find 



which proves the following statement: 



(30) 



Theorem 3 Up to local Bogolubov transformations, 
there exists a unique maximally entangled quasifree 
fermion state. 

C. Generating maximally entangled states 
Let us now consider the following one particle Hamil- 



H 











( ° 
















-In 








-1„ 








V In 












(31) 



Kp£ Here /Cab is given by the decomposition .K^ab = 
C"- © C" ® C" e C" where the first two summands are 
the one-particle space of Alice and the last two blocks 
belong to Bob. It is clear that _ff is a reflection that 
anticommutes with the complex conjugation F, i.e., we 
have H'^ = t and HT = -TH. Thus for each t G M 
the unitary operator exp(it7J) — cos(t)l 4- isin{t)H is a 
Bogolubov transformation. 

The one-particle Hamilton operator H induces a 
second-quantized Hamilton operator H that acts on the 
antisymmetrized Fock space. The relation between the 
second-quantized H and the one-particle operator H is 
given by the commutator relation 



[H,B{f)]^B{Hf) 



(32) 



which is the Heisenbcrg equation of motion. In order to 
express H in terms of Fermi field operators, we choose 
a orthonormal basis of Kab of 4n vectors e^- , e^- with 
i = 1, • • • , n such that FeJ- = e^-. The operator H can 
be written in terms of this basis according to 



>(eB+,l 



h.c. 



(33) 



The anticommutation relations Eq. (jlip can be used to 
compute commutators of the following form: 

[B(ei)tB(e2),B(/)] 

= B([F|ei)(e2|F-|e2)(ei|]/), 

where ei, 62, ./ are vectors in the one-particle space /\ab- 
This relation can directly be employed to check that the 
operator 



(34) 



H.^Y. B{e-jB{et,) + B(e+ )^-B(e-.) 



(35) 



fulfills the equation of motion Eq. ([32)) . Thus, we obtain a 
one-parameter group of Bogolubov automorphisms which 
is generated by H and which we regard now as the time 
evolution of the global fermion system. Recall, that the 
exponential form of Eq. ([5^ is given by 

exp(-iti?)B(/) exp(iti?) = B (exp(iti?)/) . (36) 

Now we apply this time evolution to the quasifree prod- 
uct state given by Eq. ((T4)) whose basis projection is 
given by Eq E. Then after time t we have created a 
quasi free state that corresponds to the basis projection 
Et = exp(— it_ff)_Bexp(it_ff). By expanding the exponen- 
tials and by making use of the relation HEqH = 1 — i?o 
we obtain 

Et ^ cos{tfEo + sm{tf{l - Eq) 

+isin(i) cos{t)[Eo,H] 

In terms of the 4x4 block decomposition, we find for the 
basis projection at time t: 



(37) 



Et := 











-iStCttn 















-iStCttn 


Cttn 





istctl„ 








S?ln 



(38) 



where ct = cos(t) and st — sin(t). Suppose Alice 
and Bob prepare their fermion systems individually such 
that the global system is in the quasifree pure product 
state with basis projection Eq. Moreover, Alice's and 
Bob's fermions interact with each other by the dynamics 
0cxp{itH)- After the interaction time t = 7r/4 the global 
fermion system is in the quasifree state that is given by 
the basis projection 





( ^n 








iln 







In 


iln 









-iln 


In 







V -iln 








In 



(39) 



which is the basis projection of a maximally entangled 
fermion state. 
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V. SPIN CHAINS AND MAXIMALLY 
ENTANGLED QUASIFREE FERMIONS: AN 
INSTRUCTIVE EXAMPLE 

In this section we consider the example of fcrmionic 
quasi-free maximahy entangled state obtained from the 
Jordan- Wigner transformation of a maximally entangled 
state in a spin chain. The spin state vector \x) defined be- 
low is related to the important class of isotropic states^i^, 
which can indeed be decomposed in terms of the identity 
and the operator |x)(xl- We will show that once written 
in terms of fermion operators, the considered state is a 
quasi-free state whose covariance matrix exactly satisfies 
the conditions stated in Theorem [2] 

We consider two blocks of spins, each block contain- 
ing n contiguous spins. By numbering the sites starting 
from the first spin of the left, the block A = {1,2,-- - ,n} 
belongs to Alice, whereas the block B = {n + l,n + 
2, • • • , 2n} belongs to Bob. We construct the state vector 
\x) as a tensor product of maximally entangled spin pairs, 
such that the {n+1 — i/)-th spin of the left block is entan- 
gled with the v-th spin of the right block, u = 1, . . . , ??. 
Fig. [T] illustrates this situation for n = 7. Explicitly, the 
vector reads 



1-iy'^n+r+iy) 



i • • • i) , (40) 



where cr^ is the usual Pauli raising operator of the j- 
th site. In the following we will drop the normalization 
factor 2-"/2. 

The Jordan- Wigner transformation is usually defined 

as 



(41) 



where Uk = cjjcfc. By substituting in the above expression 
for Ix) one gets 



\x) = n[i+^»+i-.^(^KV.]io) , 



(42) 



where s{v) ^ n„+i-^<j<„+,.(-l)"' (note that 
Cj(— = Cj). It is now possible to get rid of the 'string' 
terms s(i^) by bringing them to the right and applying 
them to the vacuum |0), so that s{iy)\Q) = |0). Indeed 
[s{iy), Cj] = whenever the site j is not contained in the 
string s(^) and this condition is always satisfied in the 
above expression for \x)- Hence all the string terms can 
be eliminated from the state, which just reads 



ix)=n(i+cn+i-.oio)=ne'"- 



'|0) . 



(43) 

The latter expression puts in evidence that |x) is a co- 
herent state (with respect to properly defined fermion 
operators), namely 



(44) 



where j^k = 1, . . . , 2n and G = iuy ® ]HI„. Here H„ is 
the trivial Hankel matrix obtained by reversing the nxn 
identity matrix, i.e., {Iln)jk = Sj+k.n+i- Note that G is 
antisymmetric and orthogonal, G^ = —G = G~^. 

The covariance matrix for a state of the form Eq. (|44p . 
where G is any real antisymmetric matrix, can be ob- 
tained from the correlations {Cjcj.) = (T — T^)jfc/4 = 
-(cjCfc) and (c]cfe) = (2 • 1 - T - r^)jfc/4, where 
T = (1-1- G) /(I — G). These expressions can be obtained 
by using the standard techniques employed to diagonal- 
ize quadratic fermion Hamiltonians (see, e.g., Ref.— ). In 
our case a simple calculation shows that T = G, so that 



(44) 



2 
2 



(45) 
(46) 



We recall that in the operators B{f) the vector / is 
ordered as f — {f^, fk^fhi fs) ^-^d that the expectation 
value of the product B{ f )B{g) is given by Eq. The 
An X An matrix S is given by 



S = 



(47) 



where the structure of the 2n x 2n submatriccs S'xy, with 
X, Y = A, B, has already been described in Sec. [Ill There- 
fore, for the Gaussian state \x) introduced above one has 



/ 1« 



(48) 



Defining J7ab = 2S'ab, it is easy to check that the matrix 
UiiB satisfies the condition FaC/ab ~ —Uiib^b- 



VI. BIPARTITE STRUCTURE AND THE 
UNTWISTING 

In this section we describe the twisted tensor prod- 
uct structure which is necessary to study entanglement 
in the global fermion algebra. The bipartite nature of 
the fermion system can be revisited in terms of the pro- 
jections Qa and Qb onto Alice's and Bob's one-particle 
spaces Kn := QaK/^b and Kb '■— Q^Khb, respectively. 
With respect to the 4x4 block decomposition, Alice's 
projection is given by 



/K 





]1„ 







(49) 



and the projection onto Bob's one particle space is 
the complementary projection Qb := 1 ^ Qa- Alice's 
local fermion algebra is generated by the opera- 
tors B{Qi^f) and, analogously, Bob's local fermion al- 
gebra ^B is generated by the operators B((5b/)- By 
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construction, the projection (and therefore Qb = 
1 — Qa) commutes with the complex conjugation which 
implies with help of Eq. (HI]): {B{QJ),B{Q^h)] = 
{TQt,f,Qiih)t = (r/,OAQB^)l = 0. Therefore, the Al- 
ice's and Bob's fermion operators mutually anticommute 
and the product of Alice's fermion operator B{Qf,f) and 
Bob's fermion operator B{Q^h) is not a tensor product, 
as it is the case for (finite dimensional) commuting ma- 
trices. 

As a matrix algebra, the global fermion algebra S'ab 
can be identified with the tensor product ® In 
the following, this isomorphism, which we denote here 
by a, is called the untwisting. It acts on fermi operators 
B{f) G 5^AB according to 

a{B{f)) := B{Q,f ) ® 1 + Ga ® B{Q^f) , (50) 

where 0a G 5^a is the reflection (selfadjoint unitary) that 
is uniquely determined by 8aB(/) — B{Rkf)Qi, and 
^A = —Qk + Qb = 1 — 2Qa is the reflection that changes 
the sign on Alice's one-fermion space Kk- We mention 
here that 0a corresponds to a string of tra-Pauli opera- 
tors, as in the case of the Jordan- Wigner transformation, 
valuating the untwisting a on a product of Alice's and 
Bob's fermion operators gives 

a{B{fk)B{h)) = B{fk)&k <E> B{h) . (51) 

Due to the occurrence of the operator Oa on the right 
hand side, the operator product B{fk)B{f^) on the left 
hand side is called a "twisted tensor product" of B{fk) 
and B{f^) which is "untwisted" by a to a tensor product 
on the right hand side. 

We note that the even part of Alice's fermion al- 
gebra is precisely the fixpoint algebra under the map 
A ^ 9k{A) := 0aA0a, i.e., an operator A belongs to 

if and only if 9k{A) = A. The analogous result is also 
true for Bob's fermion algebra. Obviously, the even part 

of the global fermion algebra consists of all operators 
that are invariant under the map F i— > QFQ = F, with 
e := OaOb. 

The essential property of the untwisting a is related 
to the odd/even nature of Alice/Bobs fermion operators. 
Namely, for each of Alice's fermion operators A € 5'a, the 
identities 

a{AB+) = A(g>B+ 

a{AB^)^ AQk^B^ . ^ ' 

are valid, where B± e is any even (odd) fermion 
operator on Bob's side. 

The twisted tensor product is hence the suitable tool to 
extend the analysis of entanglement to the global fermion 
algebra. Such an extension would be clearly desirable 
for those systems where local parity conservation cannot 
be assumed, in particular for general fermionic systems 
obtained from the Jordan- Wigner transformation of spin 
chain Hamiltonians. 



VII. CONCLUSIONS 

In this paper, we have introduced a definition of 
fermionic maximally entangled states based on the re- 
striction to the even-even part of the algebra. Such a re- 
striction is natural as far as one is interested in the phys- 
ical operators (observables) of the two considered sub- 
systems, which implies conservation of local parity. We 
have then analyzed the structure of fermionic maximally 
entangled states which are Gaussian with respect to the 
global fermion algebra. The preparation of the state is in- 
deed assumed to be constrained by global parity only, as 
for observables of the total system. Following this route, 
we have obtained a necessary and sufhcient condition on 
the covariance matrix for the Gaussian state to be maxi- 
mally entangled. We have also showed such a state to be 
unique up to local Bogolubov transformations. On the 
other hand; we have provided a unitary transformation 
which maps a Gaussian product state into a Gaussian 
maximally entangled state. An explicit example relevant 
to spin chains has also been discussed. Finally, we have 
showed how to extend the analysis of entanglement be- 
yond the restriction to the even-even part by introducing 
the twisted tensor product of fermionic subspaces. The 
latter structure is useful to study general fermionic sys- 
tems obtained from the Jordan- Wigner transformation 
of lattice spins. The investigation of this context will be 
the subject of future work. 



APPENDIX A: TRACES ON FERMION 
ALGEBRAS 

We consider the fermion algebra ^ that is generated 
by the fermion field operators B{f) with f E K = C" ® 
C". Recall that the fermion algebra can be identified 
with the algebra ^{TC) — ^ oi all linear operators on the 
antisymmetric Fock space TC = i^_(C"). On ^(H) there 
exists a unique normalized trace which is given by the 
density operator 2~"1. 

Lemma 4 The normalized trace on the fermion algebra 
^ is the quasifree state with covariance S = -2/2, i.e., the 
density operator is given by ^ 2^"_2. 



Proof: Let quasifree state that is given by the 

covariance 1/2. To prove that ^ trace, we have 

to show that the correlation functions of even products 
of fermion fields are invariant under the cyclic shift, i.e., 
we have to prove the identity 

tr(^,i/2B(/i)B(/2)---B(/2„)) 

=tr(f?i,2B(/2„)B(/i)----B(/2„-i)) 
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Let Qn C Sin be the set of all permutations q on 
{1, • • • , 2n} that fulfill for each fc = 1, 2, 3, • • • , n the con- 
straints 



g(2fc- 1) < g(2fc+ 1) 
g(2fc- 1) < g(2fc) . 



(A2) 



For our purpose, we introduce the following more handy 
representation for permutations 
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9(1) g(3) ••• q(2n-\) 



9(2) 9(4) 



(A3) 



where the first row contains the odd and the second row 
contains the even positions. The constraint for a per- 
mutation to be a member of Q„ can be reformulated in 
this representation: (1) The entries of the first row are 
increasing from left to right. (2) In each column the first 
entry is smaller than the second. By Wick's theorem, the 
left-hand side of Eq. (jAl[) can be expressed in terms of 
the two-point correlation functions by 



= 2 " ^ tq^'^]q{1k-\)-,!q(2k)) 



(A4) 



where tq is the sign of the permutation q. Let t = 
(2 3 • • • 2n 1) be the cyclic shift, then we obtain for 
the right hand side 



tr(ei/2B(/2„)B(/2)---B(/2n-i)) 



2 " X! '^9 II(r/Tg(2fc-l),/rq(2fc)) 
qeQ„ fc=l 



(A5) 



The problem is now, that the permutation Tq is not con- 
tained in Qn, namely the constraints are violated at the 
position i := q~^(2n). If we compose the cyclic shift r 
with the permutation q we obtain for Tq: 



q(l) + l ••• g(i-l) + l ••• q(2ri-l) + l 



<Z(2) + 1 



1 



q(2n) + 1 



(A6) 



The cyclic shift increases the value of the entry by one 
except the position where the value is equal to In which 
must be an even position i = q~^(2n). Thus, the entries 
of the first row are still increasing from he left to the 
right, but at position i the first entry is now smaller 
than the second. This can be cured by introducing the 
permutation <;i := [i i — 1 1 2 ••• i — 2 i + 1 ••• 2n). 
Then we build the permutation Kq := Tq(^i which can be 
represented as 



1 9(1) + 1 

<7(*-l) + l q{2) + l 



q{2n-l) + l 
q{2n) + 1 



(A7) 



By exchanging the two entries in column i/2 and shuffling 
it to the first position we have succeed in fulfilling the 
conditions (1) and (2). In fact, we have constructed a 



bijective map K:qi—f Tq<^g-i(^2n) that acts on the set of 
allowed permutations Qn- 

The product term for the permutation rq within the 
expansion Eq. (jA4p can be expressed in terms of the per- 
mutation Kq by 



fe=i 



(r/T9(2fc-l), fTq{2k)) 

n 

(r/g(,-l) + l, /l) ]^(r/Kq(2fc_l), /Kg(2fc)) (A8) 

k=2 

n 

Y\_ (r/K9(2A:-l) , fKq(2k)) ■ 



k=l 



where we have used the fact that the complex bilinear 
form (/, h) 1-^ (r/, h) = (Th, f) is symmetric. This im- 
plies the desired identity 

tiigT^^,B{f2n)B{h)---B{hn-i)) 

n 

= 2 ^ eg JJ(r/^g(2/e_l),/Kg(2A:)> 

n 

= 2-" ^^-'qY[{^fq{2k-l)Jq{2k)) (A9) 



q&Qn 



fe=l 



= 2 " ^ Eg J|(r/g(2fc-l),/q(2fc)) 

geQ„ fc=l 

= tT{gt/2Bih)B{h)---Bif2n)) ■ 

Here we have used the fact that the map k preserves the 
sign of permutations e^^ = eq. The sign of the cyclic 
shift of a set with 2n elements is e,- = —1. Namely, 
let cTfc be the transposition that exchanges k and k + 
1 we observe that r = cr2n-iO'2n-2 • • ■ ci is a product 
of 2n — 1 transpositions. Furthermore, the permutation 
Q = (ci(72 • • • (7i_i)(tTi • • • (Ti_2) can be decomposed into 
a product of 2 (i — 2) + 1 transpositions which implies 
Cc;- = —1. This yields e^q = EregEft = fg- The uniqueness 
of the normalized trace for full finite dimensional matrix 
algebras implies that the density matrix of the quasifree 
state ^'£/2 



Ql/2 = 2-"l 
which concludes the proof. 



(AlO) 



For the proof of Thcorcm[2]we need to consider all pos- 
sible traces on the even subalgebra 5^"*" which coincides 
with ^{TC+) © 3§{TC-) where TC± is the subspace corre- 
sponding to even/odd particle number. Since on each full 
matrix algebra there is - up to normalization ~ a unique 
trace, all normalized traces on are given by density 
operators of the form T = 2-"'+^{pE+ + {l-p)E^), where 
E± is the projection onto Ti±. 
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Lemma 5 Let T he the density operator for a trace on 
, then the two-point function is given by 



tr(rs(/i)B(/2)) = i(r/i,/2) 



(All) 



for all /i,/2 e K. 



Proof: Wc rewrite the density operator T asT ~ 2 " 1 — 
2-"+i(2p - 1)6 where 9 = E+ - E_. Let Ci,c| be 
the standard creation and annihilation operators on Fock 
space and let TV = c^Cj be the particle number op- 
erator. There exists an orthonormal basis (e/) of eigen- 
vectors, i.e., Nei = njej where there are (^) eigenval- 
ues with ni = k. To be more explicit, we choose a or- 
thonormal basis for the one-particle space (e,;|i = 1, • • • n) 
and associate to each subset / C N := {!,■■■ ,n} 
we associate the vector ej := e^^ A e,;^ A • • • e,;^. with 
/ = {ii, • • • ^i]^} and ii < i2 < ■■■ < ik- Here 
i/j /\(f) :~ (g) (/) — (g) V') is the antisymmetric tensor 

product. The standard creation operator acts according 
to 



cUi ■■= A e.i^ A e,;^ A • ■ 



(A12) 



which vanishes if j El. li j G I and let p{j\I) be the 
number of elements in / that are smaller than j, then we 
get 



cUi = {-ir^^^'^d[jeN\I]eiuj 



(A13) 



where we use the convention (5 [true] ~ 1 and (5[f alse] 
0. Obviously, we get for the annihilation operator: 



Cjei = {^1)pW\=)S[j e I]ej\, 



(A14) 



Indeed, (e/) is an orthonormal basis of eigenvectors 
with Nej = \I\ej. Since 6 = (-1)^, it follows that 
{ei,QciCjei) ~ (— l)l^l(e/, CiCje/) ~ 0. For the term 
{ei,Qclcjei) we find 



Thus the trace of Ocjcj can be computed according to 
tr(eclc,)= ^(-l)l^l5[ie/]<5,, 

JCN\i 

Here we have used the fact that N \ i contains as many 
even-elementary subsets as odd-elementary subsets, with 
an even. This implies ti{QB{fi)B{f2)) = and with the 
help of Lemma [4] we find the desired result: 



tr(rB(/i)B(/2)) = 2-"tr(B(/i)B(/2)) = -(r/1,/2) 



(A17) 



With help of the two lemmata above we prove the fol- 
lowing theorem: 

Theorem 6 There is a unique quasifree normalized trace 
on the even subalgebra which is determined by the 
covariance i/2. 



Proof: By Lemma IH the quasifree state with covariance 
1/2 is the normalized trace on the full fermion algebra 
g'. The restriction to the even subalgebra 3"+ is again 
a quasifree normalized trace. On the other hand, each 
state T on ^ that vanishes on the odd part ^~ and 
which induces a trace on the even part g'^ has the form 
T = 2-"l + 2-"+i(2p - 1)9. If T is quasifree, then it 
is uniquely determind by the two-point correlation func- 
tion. Suppose T = gs is quasifree with covariance S, 
then it follows by Lemma [5] that 8—1/2. ■ 



(e/,ectc,e/) = (-l)!^l<5[ze/]4 



(A15) 
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